% documentation for grun.data
% (Grunfeld investment data)

% variables I,F,C for firms:

% general electric
% westinghouse
% general motors
% chrysler
% us steel

% order is i,f,c (see below for more details)

clear all;

load grun.dat;  % grunfeld investment data
                % see page 650, Green 1997

y1  = grun(:,1); %   I general electric
x11 = grun(:,2); %   F
x12 = grun(:,3); %   C

y2  = grun(:,4); %   I westinghouse
x21 = grun(:,5); %   F
x22 = grun(:,6); %   C

y3  = grun(:,7); %   I general motors
x31 = grun(:,8); %   F
x32 = grun(:,9); %   C

y4  = grun(:,10); %   I chrysler
x41 = grun(:,11); %   F
x42 = grun(:,12); %   C

y5  = grun(:,13); %   I us steel
x51 = grun(:,14); %   F
x52 = grun(:,15); %   C

% set up a structure for y containing y's for each eqn
% (order follows that in Green, 1997)
y(1).eq = y3; % gm
y(2).eq = y4; % chrysler
y(3).eq = y1; % general electric
y(4).eq = y2; % westinghouse
y(5).eq = y5; % us usteel

n = length(y1);
iota = ones(n,1);

% set up a structure for X in each eqn
X(1).eq = [iota x31 x32];
X(2).eq = [iota x41 x42];
X(3).eq = [iota x11 x12];
X(4).eq = [iota x21 x22];
X(5).eq = [iota x51 x52];

vname1 = ['I gen motors  ',
          'constant      ',
          'fgm           '
          'cgm           '];
          
vname2 = ['I chrysler    ',
          'constant      ',
          'fcry          '
          'ccry          '];

vname3 = ['I gen electric',
          'constant      ',          
          'fge           ',
          'cge           '];
          
vname4 = ['I westinghouse',
          'constant      ',
          'fwest         '
          'cwest         '];
          
vname5 = ['I us steel    ',
          'constant      ',
          'fuss          '
          'cuss          '];


% do sur regression with iteration
neqs = 5;
iflag = 1;
result = sur(neqs,y,X,iflag);

% no iteration done here
result2 = sur(neqs,y,X);


vname = [vname1
         vname2
         vname3
         vname4
         vname5];

prt(result,vname);
prt(result2,vname);